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A DECOMPOSITION FOR THE SCHRODINGER EQUATION WITH 
APPLICATIONS TO BILINEAR AND MULTILINEAR ESTIMATES 

FELIPE HERNANDEZ 


Abstract. A new decomposition for frequency-localized solutions to the Schrodinger equa¬ 
tion is given which describes the evolution of the wavefunction using a weighted sum of 
Lipschitz tubes. As an application of this decomposition, we provide a new proof of the 
bilinear Strichartz estimate as well as the multilinear restriction theorem for the paraboloid. 


This paper introduces new way of decomposing solutions to the Schrodinger equation, 
which approximates the evolution of the mass distribution of a solution as a sum of space- 
time Lipschitz tubes. There are several ways in which this decomposition differs from existing 
decompositions of solutions to the Schrodinger equation, such as the wavepacket decompo¬ 
sition. First, whereas the wavepacket decomposition breaks the solution up into tubes of 
width i? 1//2 and length R, the tubes appearing in this decomposition have unit width and 
can be made arbitrarily long. Second, the tubes involved are curved instead of straight. 
Finally, the decomposition involves no cancellation. Whereas in the usage of the wavepacket 
decomposition it is often necessary to induct on scales in order to resolve possible cancella¬ 
tions between tubes, all pieces of the decomposition introduced here are positive so there is 
no opportunity for cancellation. The lack of cancellation to exploit and the relatively loose 
control of the shape of the tubes means that this decomposition seems unable to prove many 
dispersive inequalities used in the literature on the Schrodinger equation. However, the de¬ 
composition is useful enough to prove the bilinear Strichartz estimate and a special case of 
the multilinear restriction estimate. One motivation for this work was to find proofs of these 
estimates that can be generalized to non-Euclidean settings such as hyperbolic space. 

Now let us set up the general situation for our results. Let u : x M —* C be a solution 

to the Schrodinger equation 

(1) id t u + A u = 0 

with initial condition u(x, 0) = Uq(x) G L 2 (R d ). For convenience we will write u t for the 
function Ut(x) = u(t,x). If u t is localized in frequency space, so that the Fourier transform 
supp u t C Bi is contained in the unit ball, then physical intuition about the Schrodinger 
equation suggests that u t has two important properties. First, because of the uncertainty 
principle, one expects that u t (and thus \u t \ 2 ) does not vary much (is “locally constant”) 
on unit scales. Moreover, because u t has bounded momentum, one expects that u t enjoys 

Date: July 28, 2015. 

2010 Mathematics Subject Classification. Primary 35Q41; Secondary 42B37. 

Key words and phrases. Schrodinger equation, bilinear Strichartz, multilinear restriction. 

1 



2 


F. HERNANDEZ 


something akin to finite speed of propagation. Together these suggest that it is possible 
to describe the evolution of the mass \u t \ 2 in terms of discrete “packets” that travel along 
finite-speed paths. The point of this paper is to make this intuition precise. 

To set up some more notation, suppose 7 : M —> is some differentiable path which will 
represent the motion of one such packet. We say that 7 has speed at most V if | V7I < V. 
Moreover we define the tube of width r centered at 7 by 

T 1>r := {(ay t) G R d x M; |x — 7(f)| < r}. 

We will use T 7jr and its indicator function interchangeably. Then our main result is the 
following decomposition. 

Theorem 1 (Skinny Lipschitz Tube decomposition). Let u t solve the Schrodinger equa¬ 
tion (HD with frequency localization suppwo C B\. There exists a radius r, a speed limit 
V, and a constant C with the following property: For any R > 0 , there exists a countable 
collection of paths {7*} with speed at most V and weights {wi} such that 

(2) \u t (x )\ 2 < Yw iTl Ax,t) 

i 

for any (ay t) e R d x [— R, i?], and moreover this bound is efficient in the sense that 

( 3 ) Y r d Wj <c[ \u t \ 2 dx. 

jR d 


Remark: By rescaling the weights tty so that JT-uy = 1 > it is possible to think of them 
as probabilities. In this interpretation, wy is the probability that the particle represented by 
the wavefunction ut takes the path described by 7. However, since the paths are far from 
unique, this is not a very accurate description. 

The proof of this theorem has two main ingredients. The first ingredient is a pair of 
inequalities which demonstrate two key properties of u t . The first inequality says that u t is 
locally constant , in the sense that \u t (x )\ 2 can be bounded pointwise in terms of a weighted 
average centered at x. The second says that \u t \ 2 has finite speed. The statement of the finite 
speed property is slightly more complex, but suffice it to say that this is where we encode 
the “bounded momentum” intuition. The second ingredient of the proof is combinatorial. 
We state and prove an analogue of Theorem (T] in a discrete situataion. The main idea here 
is to construct discrete paths one step at a time. The existence of these one-step paths is 
guaranteed by the Max-Flow Min-Cut theorem. 

By itself, Theorem |Tj is unsuitable for applications which require a variety of frequency 
localizations. This is easily fixed with an application of scaling and Galilean symmetries. 

Corollary 1 . Let ut be a solution to the Schrodinger equation with frequency localization 

suppfi 0 C B p (g), 


A DECOMPOSITION FOR THE SCHRODINGER EQUATION 


3 


where f E R d is some average momentum, and p is the frequency uncertainty. Let r,V,C> 0 
be the same constants given in Theorem |7J and let T > 0 be arbitrary. Then there exists a 
set of paths with (Vy* — £| < Vp, and weights such that 

hwi 2 < E (x,t) 

i 

and 

(4) 'V](rp~ 1 ) d w i <C [ \u t (x)\ 2 dx. 


Proof of Corollary [7j using Theorem [/} Consider the function 

u'(x,t ) = e l ^ 2t / p2 - £ -' x / p ' > u(xp- 1 + tfp~ 2 ,tp~ 2 ), 

which is still a solution to the Schrodinger equation (by Galilean and scale invariance), and 
which moreover satisfies the frequency localization supp u' C B\. Apply Theorem [1] to u! 
and then apply the scaling and Galilean symmetries again in reverse. □ 

As mentioned above, Theorem [Q and Corollary [Tj seem unsuitable for some applications. 
In particular, one apparent obstacle is that Theorem |T] provides no guarantee about the 
dispersive behavior of the tubes. It may be that there is a single tube which remains coherent 
for all time. Thus, Theorem |Tj cannot by itself prove something such as the Strichartz 
inequality. However, there are examples of useful estimates in which dispersion does not 
play as large a role. The two we illustrate here are the bilinear Strichartz estimate and the 
multilinear restriction theorem. 

0.1. Bilinear Strichartz. The bilinear Strichartz estimate was introduced by Bourgain 
in [ Bou98j . It handles the interaction between high-frequency and low-frequency solutions 
to the Schrodinger equation. To state it, let 

Av := {£; N/2 < |£| < 2N} 

denote the annulus at scale N. 

Theorem 2 (Bilinear Strichartz [Bou98] ). Let M N, and let u 0 ,v 0 E L 2 (R d ) have the 
frequency localizations suppu C A N , supp v C A M . Then 

M ^ d ~ 1 )/ 2 

(5) ll^ilUl^IRdxIR) < C jy 1/2 —11 ^0 11 (jjd) 11U 0 11 L| (R d ) • 


This estimate is one of the ingredients in the I-method, which has been used to prove 
global well-posedness for the nonlinear Schrodinger equation for rough data |CKS + 02j . The 
first proof of the bilinear Strichartz estimate relied heavily on properties of the Fourier 
transform which do not easily generalize to non-Euclidean settings. Since the first proof by 
Bourgain, Tao has given a more physical-space proof in [TaolOj . Moreover it seems that 
the wavepacket approach to the bilinear Strichartz estimate for the wave equation described 
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in [ KRT02 ] might work for the Schrodinger equation as well. But to my knowledge, the 
only proof of the bilinear Strichartz estimate on a manifold is due to Hani [Hanl2a| . Hani 
used this estimate and the 1-method to prove a global well-posedness estimate for the cubic 
nonlinear Schrodinger equation in [Hanl2bj . The approach in this paper to the bilinear 
Strichartz estimate is quite different from any mentioned above, and I hope that it may be 
extended to non-Euclidean settings such as hyperbolic space. 

0.2. Multilinear Restriction. The multilinear restriction estimate was first stated by Ben¬ 
nett, Carbery, and Tao [BCT06] . It is a version of the restriction problem which is simplified 
to handle only quantitatively transverse interactions. Here we work with a very special case, 
since we only want to show how the decomposition of Theorem [Tj can be used. 

Theorem 3 (Multilinear restriction |BCT06j ). Let Ui be solutions to the Schrodinger equa¬ 
tion for 0 < i < Let r = (lOcfR 1 and suppose that uq has frequency support suppwo C 
B r (0) and the Ui have frequency support supp u t C Brief), where e, is any orthonormal basis 
for W l . Then for every e > 0 there is a constant C £ > 0 such that 

d d 

<R i= 0 i=0 

for any R > 0. 

This kind of estimate has been used to make progress on the restriction problem jBGll] 
and more recently to prove the £ 2 decoupling conjecture [BD14j . We remark that the main 
difficulty in proving the multilinear restriction is an incidence geometry problem about tubes 
(see [Gutl5j for a short discussion of the problem), and our only contribution is to provide 
a new derivation of the restriction estimate from the geometric problem. In particular this 
derivation avoids the use of induction on scales. 

0.3. Plan of the Paper. The paper is organized as follows. Section Q] contains a derivation 
of the locally constant and finite speed properties to frequency localized solutions. Section [2] 
proves a discrete analogue of the decomposition. In Section [3j the full proof of Theorem |Tj 
is given. Section [4] proves the bilinear Strichartz estimate, Theorem [2l and the multilinear 
restriction estimate, Theorem [3l using the tube decomposition in Corollary [T] 

Acknowledgements. I would like to thank Larry Guth for useful comments. I also thank 
Vincent Tjeng for showing me how to use the Max-Flow Min-Cut theorem to prove Lemma[4j 

1. Properties of Frequency-Localized solutions 

In this section we prove variants of the locally constant (LG) and finite speed (FS) prop¬ 
erties for solutions u t to the Schrodinger equation satisfying 

(7) supp Ug C B|. 

Tn the statement of this theorem we drop the notation that Ut.(x) = u{t,x). 


/ / 

J |x|<tt J id 
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1.1. Locally Constant. We begin with the locally constant property. We say that a 
bounded function / G L°°(M d ) is localized to the unit ball if there exists some constant 
C such that \f(x)\ < C and 

|/(x)| < C\x\~ 10d , |V/(x)| < C|x| _1M_1 . 

Proposition 1. There exists some positive function fi which is localized to the unit ball and 
such that for any function u 0 satsfying the frequency j7J) , 

(LC) sup \u 0 (y)\ 2 < C [ \u 0 (y)\ 2 p(x - y) dy. 

yGBi(x) J 


Moreover, // can be chosen such that the translates of ji form a partition of unity, 

(8) M x — a) = !, 

a£Z d 


and such that 

(9) c\x\- 10d < p. 


Remark: The extra conditions (J8J) and (jUJ) can be tacked on by starting with some p" 
localized to the unit ball and then making it larger. At first glance this seems quite lossy; 
the point however is that in the proof we will use the bound (ILCD only once, and the extra 
mass in p will allow us to arrive at a finite speed property (in the next subsection) which 
has no loss of constants, and which will be iterated many times. 

Remark: The bound in (J9]) could be tweaked in several ways. For the purposes of the next 
section, it would suffice to know that /i is smooth (approximately locally constant on unit 
scales) and that 



\dM < C 



R 


whenever E cR d satisfies EnBi =0. This is a more general condition which may be useful 
in non-Euclidean settings, but the easiest way I know to verify it is to enforce that p has 
polynomial decay, hence the bound (J9j). 


Proof. Let y G L 2 (M d ) satisfy y(£) = 1 for \f\ < 1 and y(£) = 0 outside |£| > 2, with a 
smooth cutoff in between. By choosing the function appropriately, we can enforce that y 
localized to the unit ball. The Fourier localization of Uq yields Uq = y * m 0 . This allows us 
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to bound, using Cauchy-Schwartz, 


\u(x)\ 2 = \x*u(x)\ 2 = 


u(y)x(x-y) dy 


< / \x(x-y)\dy) / \u(y)\ 2 \x{x - y)\dy 


= c d \u{y)\ 2 \x{x-y)\dy. 


Now define the function /i" by 

y"(y)= sup \x(z)\- 

\z-y\<x 

The fact that x is localized to the unit ball implies that y” is also, with a larger constant. 
We now modify y" to satisfy (J8J) and (J9]) . As long as these modifications only increase y',we 
can ensure that (1LCD is still satisfied. 

First, since y" is localized to the unit ball, we can find some function y! which is smooth, 
is localized to the unit ball, has polynomial decay y' > c|x| _10d , and satisfies y" < y'. 
Moreover, each dilation y' R (x) = y'fx/R) is still localized to the unit ball and satisfies (JH]) 
with different constants. By making R large enough the periodic function 


P ( x ) =^y(x- a) 

aez d 

becomes very smooth, so that y = y'(x)/p(x) which satisfies ([8j) still obeys all previous 
bounds. □ 


1.2. Finite Speed. In this section we prove a version of the finite speed property for fre¬ 
quency localized solutions to the Schrodinger equation. We first establish some notation for 
the section. Define the set 

(10) H := {h £ Z d ; max \hf\ < 1}, 

l<i<d 

which has cardinality 3 d and is the cube centered at the origin. Given any set A C Z d , we 
let A + H denote the sumset 

A + H := {a + h G Z d ; a G A, h G H}. 

In other words, A + H consists of all lattice points in Z d which are within a distance 1 of A 
in the £°° norm. We also define dA = (A + H) \ A, the lattice points that are within 1 of A 
but are not in A. 

We will also be using y and its translates from Proposition [Q For convenience, let y a (x) = 
y(x — a ) be the function y centered at the lattice point a £ 7L d . Moreover, for any set B C Z d , 
we write y B = Y,beB Pb- 

Proposition 2. Let uq have frequency localization as in © and u t = e ltA uo be the free 
Schrodinger evolution of u 0 . Moreover let y satisfy the conclusions of Proposition^ 
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Then there exists some time r > 0 and such that for any A C Z d and 0 < t < r, 


(FS) 


\u t \ 2 HA < J \Uo\ 2 hA+H, 

KIVa < / \u t \ 2 Ta+h- 


Remark: By time reversal symmetry, the equations in (IFSI) are equivalent, so we only 
focus on proving the first. 

We will need a few lemmas in order to complete the proof of this proposition. 

Lemma 1. There exists a constant C > 0 such that for any A cZ d and any 1 < j < d, 
(11) \djix A \ <C\^ dA \. 


Proof. Let x G M d , and let dist(x, dA) denote the minimum distance between x and a point 
in dA. We split the analysis into cases. The first case occurs when x is near the boundary 
of A, so that dist(x,cM) < lOOd. Then /j,qa(x) = f^(l), and \dj/j,A\ = Od( 1), so the bound 
holds for some constant. 

The second case is when x is in the “exterior” of A, so dist(x, A ) > lOOd. Then we write, 
expanding the definition of jua, applying the triangle inequality, and using the fact that /! is 
localized to the unit ball 

\djHA\ < C \x — a | -10d-1 . 

aeA 

We would like to turn this sum into an integral. Let Q(A) C denote the union of unit 
cubes centered at the points of A. Since we are well away from the origin, the function 
| 2 /|—i ° rf —1 doesn’t change much on unit scales, so 




| — 10 < 2—1 


< c 


\x — a 


- 10 ( 2-1 


da. 


aeA 


IQ(A) 


Now observe that the integrand can be written as — 9ddi(xi\x\ wd x ), so by integrating by 
parts we obtain 


\x 


— a| _10d_1 da < c 


\x — a 


10(2 


Iq(a) 


' dQ(A) 


dcr(a), 


where da denotes the surface measure on the exposed faces of the cubes. We can assign each 
face of an exposed cube to a point in dA, with each point in dA being chosen at most 2 d 
times. Again because we are away from the origin, the function |i/|~ 10d doesn’t change much 
on unit scales, 


\x — a 


- 10(2 


'dQ(A) 


da (a) < C E \x — 


1 - 10(2 
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Applying m and chaining the inequalities we can conclude that 

Sj pa A C* I I • 

Finally we must handle the case that x is in the “interior” of A, so that dist(x, A) < lOOd. 
This time the bound 

\djHA\ < C \x — a | _10d_1 

aeA 

is much too lossy because of the cancellation involved. To exploit this cancellation, we 
observe that because p forms a partition of unity, 

I hA | \9jpA c |; 

where A c = 7h d \ A is the complement of A. Now dist(a;, A c ) > lOOd, so apply the same 
argument as above. □ 

Lemma 2. Let K £ S( R d ) be a Schwartz function (smooth and rapidly decaying). Then 
there exists C > 0, independent of K, such that for any B C 7L d , 

(12) \K * n B (x)\ < C (^J \K(y)\(l + \y\ 10d ) dy^j p B (x) 


Proof. By the triangle inequality and translation invariance, it suffices to show that 

|K * p(x)\ < Cp(x). 

Now let y £ W d . We claim that 

(13) p(x-y) <C(l + \y\ 10d )p{x). 

Again we split the analysis into cases. The first case occurs when |a:| < 50|r/| and \y\ > 10. 
Then p(x) > c\^t)y\~ wd 1 so 

y(x-y) <C < C\y\ 10d p(x). 

The second case is |x| < 50|?/| and \y\ < 10. Then both p(x) and p(x — y) are 0(1). Finally, 
if |x| > 50|?/| then 

y(x - y) ~ p(x) 

because both are 0(|x| _lod ). This concludes the proof of the claim (fT3lh 
We use this claim to bound the convolution: 

I K * p(x)\ = | J K{y)p{x - y) dy\ 

< J \K(y)\\n{x ~y)\dy 
<Cp(x) J \K(y)\(l + \y\ 10d ) dy. 


Now because K is rapidly decaying, the integral on the right is bounded by some constant 
so we are done. □ 
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Proof of Proposition® Subtract \u 0 \ 2 pa from both sides of (IFSD . apply the fundamental 
theorem of calculus, and use the local conservation of mass 

d t \u t \ 2 = dj Im (u t djU t ) 

to note that it suffices to show 


d s \u s \ 2 ds p A = 


dj lm(u s dfu s ) ds) p A < / \u Q \ 2 p dA - 


(\u t \ 2 -\u 0 \ 2 )p A = 

J J Vo / J Vo 

We may integrate the LHS by parts. Moreover we use u s = e lsA u = e tsA (x*u) = (e lsA x) *u, 
and define K s = e lsA x■ We can rewrite this inequality as 


Im( 



(K s * uo)(djK s * uo)djp A ) < / \u 0 \ po a- 


We now need a few facts about K s and dK s . They are both concentrated in the ball 
|x| < s and have rapidly decaying tails (the estimate on the tails is a repeated application 
of integration by parts). In particular, the quantities 


J \K a (ym+\yr)dy, J \djK s (y)\(l + \y\ Wd ) dy 

are uniformly bounded near s = 0. We can therefore apply Cauchy-Schwartz several times 
to bound the LHS by 



Im( / / (K s * u 0 )(djK s * u 0 )djp ) < (\K S \ * |w 0 |)(|<V’C s | * \u 0 \)\djfi A \ 



< C 


\ 1/2 

(|/i s | * \uv\Y\djPA\ 




i\djK s \ * \u 0 \) 2 \djp A \ 


< Ct sup { / |n 0 | 2 |^| * \djP A \, / |n 0 | 2 |^iL s | * \djp A \}. 
0<s<T,l<j<d J J 

By applying Lemma [T] followed by Lemma ® we can bound the right hand side by 

Ct [ |-Uo| 2 |/pm|- 


Upon taking r < 1/C we are done. 


□ 


2. The Discrete Situation 

In the continuous situation, we have an evolving mass distribution \ut \ 2 and we would like 
to understand it in terms of packets moving around according to Lipschitz paths. This can 
be discretized in space and time to come up with the following set-up. We have a (possibly 
infinite) digraph G = (V, E) with bounded degree (self-loops are allowed). The nodes 
represent a discretization of the space (in our case V = Z d ), and edges represent possible 

2 In our notation, a digraph has directed edges but it is allowed for both the edges (u,v) and (v,u) to 
exist. 


1/2 
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movements of mass within a timestep (in our case (a, b) G Z 2 is an edge if ||a — 6 ||< 1). 
Instead of an evolving mass distribution \u t \ 2 there is a sequence of mass distributions w\ : 
V —> R + for i — 1, 2, • • • , N. We say a sequence of vertices p = (p(l), • • • ,p(N )) is a path 
of length N if {p[i),p[i + 1)) is an edge for each 1 < i < N. The set of all such paths is 
denoted Vn■ We want to describe this evolution of mass as a sum of packets that each move 
according to different paths. 

Question 1. Is it possible to find a weighting a : Vn —> R + on the set of paths Vn of length 
N such that for all 1 < i < N, and all vertices v G V, 

wfiv) = Y “(p) ? 

p&V;p(i)=v 

Sometimes it is impossible to come up with such a weighting. For a trivial example, 
there may be no edges in E at all, and so there can be no paths. The only allowable mass 
distribution is the trivial one, w t = 0. A less trivial constraint is that Wi must satisfy a 
conservation of mass. Indeed, observe that 

Y Wi a (p) = 5Z a (p) 

vE.V vE.V p€'P;p(i)=v p(zP 

is independent of i. In fact, Wi must satisfy a local conservation law. To state it, we make a 
few definitions. For A C V define A^ + (A), the outgoing neighborhood of A, to be the set 

N + (A) := {v G V] ( a,v ) G E for some a G A}. 

Similarly define the incoming neighborhood, 

N~(A) := {vGl'; (v, a) G E for some a G A}. 

By an abuse of notation we will write A^ ± (u) instead of A^ ± ({u}) when A consists of a single 
vertex v. 

The key property here is that if p is a path and p[i) G A, then p{i + 1) G iV + (A) (and 
p(i — 1) G N~{A)). Thus, 

Y w ^ = 5Z a (p) - a (p)= Wi ( 6 )- 

aeA p(i)eA p(i+l)£N+(A) b£N+(A ) 

Likewise _ __ 

Y v h(a) < Y 

aeA beN-(A) 

The answer to our question above is that these local conservation laws are not only necessary 

but also sufficient. 

Proposition 3. Let G = (V, E) as above, and let Wi : V -A- M + be a sequence of positive 
weights on V. Suppose that for every A C V, and 1 < i < N the local conservation law 
holds, 

(14) J^(a)< Y w i+ i (b), 

aeA b€N+(A) 
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and in addition the global conservation law holds, 

( 15 ) '%2w 1 (v) = '^2w n (v). 

vGV vdV 

Then there exists a weighting a : Vn —> M + on the set of paths of length N such that for 
every 1 < i < N and every v G V, 

( 16 ) Wi(v) = a (P )■ 

p(i)=t) 


First we prove the Proposition in the case N = 2. Then we deduce the case of general N 
from the N = 2 case. 

2.1. The case N = 2. When N — 2, Proposition [3] is a consequence of the Max-Flow 
Min-Cut (MFMC) theorem. An introduction to the MFMC theorem can be found online, 
for example see [ Jos07 . ISM) . The theorem is originally due to Ford and Fulkerson [FF56j . 
Usually this theorem is stated in the case of a finite network, but here our vertex set V 
is infinite. In general this can be problematic, but the weights we consider are absolutely 
summable so there is no real difficulty. 

Now we set up notation so that we may state the MFMC theorem. Let U be a set of nodes 
with two distinguished elements, s and t, which represent the source and sink respectively. 
Let c : U x U —* M + . We say c is a capacity function if c(u,v ) > 0 implies c(v,u ) = 0. That 
is, if we think of c as being the capacity of a network of pipes connecting the nodes of U, 
the pipes must only go in one direction. A cut is a subset S C U such that s G S and t fL S. 
Given a capacity function c, we define the capacity of the cut S, written cap S, by 

cap S'= c(u,v). 

u£S,v<£S 

The dual object of a cut is a flow. A flow is a function / 

Kirchoff’s laws 

fin(v) := = £/(„,«) =: 

u€U u€U 

for all v G U \ {s, t}. The value of a flow is given by 

val / — £/M. 

u&U 

Because / satisfies Kirchoff’s laws except on s and t, we may also write 

val / = £/(«,*). 
ue u 

The MFMC theorem relates the maximum value of a flow to the minimum capacity of a cut. 


: U x U —> M + satisfying the 

fout ( t ) 
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Theorem 4 (Max-Flow Min-Cut |FF56j ). Let U be a finite vertex set with source s € U 
and sink t e U, and let c : U x U —> M + be a capacity on U. Then there exists a flow 
f : U x U —>■ M + with f < c such that 

val f = min cap S. 

scu y ’ 

where the minimum ranges over cuts S C U with s € S and t ^ S. 


Notice that flows, which are weights on edges, can also be thought of weights on the set of 
paths of length 2. This provides the main connection between MFMC and the N = 2 case 
of Proposition [3] To illustrate this connection better, we set forth some definitions. Let V\ 
and V 2 be two identical copies of V, and let U = {s, i} U V\ U V 2 . Let c : U x U —> M + be 
the capacity function defined by 


(17) 


Wi(vf) 

if u — s, v! e Vi 

Wi(u) 

if u E Di, u! e N + {u 

w 2 (u) 

if u e V 2 , v! — t 

0 

else. 


If we could apply MFMC directly to U with capacity function c, then the local conservation 
laws GW) and (fT51) would ensure that there is a suitable flow. 

Unfortunately U as described may be infinite. We therefore need to apply MFMC to finite 
subnetworks of U. Let A C V be a finite subset of V, and let U a = {s, t}UT 1 UA^ + (A 1 ) C U, 
where we think of A\ c V\ and A^ + (A!) C V 2 . We can also define the capacity ca = c\u A , 
which is simply the restriction of c to Ua■ We can apply MFMC to the network described 
by the nodes Ua and capacity ca to obtain the following result. 


Lemma 3. Let A C V, and let Ua and ca be as described above. Then there exists a flow 
/a : Ua x Ua —> such that /a < ca and 

/a(s,u) = wflu) 


for every u £ A. 


Proof. Consider the flow network on Ua with capacity ca- Let S C Ua be a cut. We can 
write S = {s} U Si U S 2 , where Si C Ai and S 2 C N + (A). Now we write cap S using the 
definition (jTTJ) of c 


capS= c(s,u) + c(u, v ) + c(v, t ) 

mGAi\Si mGSi V&S 2 


uGS 1 

veN+(Si)\Si 


w i ( m )+ wflu) + ^2 w 2 (v). 

u(E.Ai\Si uESi v£S2 


1 

ceN+(S 1 )\ 5 2 


Define a subset S[ C Si by 


S[ := {u G Si, N + (u) C S2}. 
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If u G Si\ S[, then there is some v G N + (u ) such that v (f S 2 , so that c(u,v) = Wi(u) 
contributes w\(u) to the sum in PIS'. Moreover we have S 2 C N + (S[), so 

cap S'> Wi(u) + w 2 (v) 

ugA 1 \S[ vGN+(S[) 

> Y w i(u) + Y w i( u ) = ^2 w i( u ) 

ugAi\S[ uGS[ ugA 1 

where we have used the local conservation law (fT4l) with i — 1. The expression on the RHS 

is exactly cap({s}), so by MFMC there exists a flow f a on Ua with f a < ca such that 

val Ja = Y Wl ( u )- 

ugA 

On the other hand since f a(s,u ) < c(s,u ) = Wi(u), 

val /a = YfA( s > u ) - '22 w ^( u )- 

u£A u£A 

Thus / a(s,u ) = u>i (u) for all u e A, as desired. □ 

A compactness and limiting argument allows us to make a conclusion about the entire 
network U. Formulated in slightly different notation, we obtain the N = 2 case of Proposi¬ 
tion [3] 

Lemma 4. Let G = (V, E) be as in Proposition [U and let Wi and vj 2 be two weights on V 
satisfying (UTD with i = 1 and (1T31) for N = 2. Then there exists a function f : E —$■ M + on 
the edges such that 

wi{v)= Y /(o«) 

uGN+(v) 

and 

w 2 (v) = Y f( U ’ V )' 

u£N— ( v ) 


Proof. Let A k C V be an increasing sequence of sets, Ai C A 2 C • • • , such that for every 
v G V there is some k such that v G A k . Let Ua be the network described above. For each 
k, Lemma [3] provides some flow f k such that f k < c and f k (s,u ) = W\[u ) when u G V. 
We show that, up to subsequence, f k converges to a flow / on Ua in U. This is done by a 
compactness argument. 

Indeed, let e > 0. We will prove the existence of a subsequence k( 1), k( 2), • • • such that 

\fk(i) — fk(j) | < £ 

(u,u')GUxU 


( 18 ) 
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for all i, j > 0. First choose K so large that 

Y Wl ( M ) - 

uG V\A k 

which is possible since W\ is absolutely summable and A is an increasing sequence which 
converges to V. 

There are finitely many pairs (u,u') G Ua k x Ua k , so by compactness we can find a 
sequence k( 1 ), k(2 ), • • •) such that 

(19) Y \fk(i)(u,u') - f Hj) (u,u')\ < £ 

(u,u')eU a k x Ua k 

for all i,j. Now we write 

^ ^ I ("U, M ) | y ^ | fk(j) (®, U) \ 

(u,u')eUxU u^Ak 

+ Y ~ fk(i)(u,v)\ 

u£A k 

vGV 2 

+ Y \fk(i)( v it)~ fk{i)(v,t)\ 

v£N+(A k ) 

+ Y I fk(i)(u, v ) - fk{j)(u, v )\. 

uGAk 

vGV 2 

We bound each term separately. Because of (1191) the last term is bounded by e. The first 
term is bounded by 2e due to the triangle inequality, the fact that / < c, and the definition 
of Ax- The next two terms are actually identical to the first, once you apply the triangle 
inequality and the Kirchoff’s laws for f^u) and fk{j)- We have therefore shown (fTHD . 

A standard diagonalization argument now shows the existence of a subsequence, which we 
will not relabel explicitly, which is Cauchy in E 1 . More precisely, this means for every e > 0 
there exists M such that when n, m > M , 

Y I fn( u , u ') ~ fm{u,u')\ < £. 

(u,u')GUxU 


By the completeness of E 1 , this converges in E 1 to some function / : U x U -A- M + . This 
convergence is strong enough to ensure that / is still a flow, and moreover f(s,u ) = W\{u) 
for all u G V. This, combined with the Kirchoff’s laws, implies that 

wi (u)= Y /( M W) 

vGN+(u) 
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for all u G V. Moreover, Kirchoff’s laws and the global conservation law (fT51) ensure 

^ w i(«) = val f = Y = Y Wl ^' 

uev v£V2 vev ugv 

Since f(v,t ) < w 2 (v) pointwise, this equality enforces f(v,t ) = W 2 (v) for all v G V. We have 
therefore shown that / has the desired properties. 

□ 

2.2. The case N > 2. This short section just shows how to iterate the N = 2 case to 
complete the proof of Proposition [3J 

Proof of Proposition^ We start with a simple observation. By (fl4jl . since N + (V) C V, 

Y^ w i( v ) < < • • • < w N (v). 

vev vgv vgv 

But by (TT5j) all of these inequalities are actually equalities. 

Now we consider the graph on the vertices Vjjy] ■— V x [N], where (u, i ) G V[jv] is connected 
to (v,i + 1) G V[n] if v G N + (u). Each edge can be uniquely assigned a triple ( i,u,v) 
where 1 < i < N and v G N + (u). Denote by Epv-i] the set of such triples. We write 
V[n] = W 2 U-■ -U Vn where Vj = V x {f}. Similarly we write E[ N _ i] = EiUE 2 U- ■ -UE n _i, 
where each E t has the first coordinate of its triple equal to i. We would like to assign a 
weighting function / : Ey tv— i] —> 1R + to each of these edges such that for each 1 < i < N and 
veV, 

(20) Wi(u)= /(b M w), 

veN+(u) 

and for each 1 < i < N, 

( 21 ) Wi(v) = Y /(*“ !,«,«)■ 

u£N~(v) 

This can be done by finding, for each 1 < j < N, a flow fj : E t — y M + such that f! 20 p is 
satisfied for i—j and (|2T1) is satisfied for i — j + 1 . But the existence of fj is exactly what 
is guaranteed by Lemma 01 Now simply define / on Ey tv— i] by f\ E = fj- 

Finally we use / to define a weight on each path. We approach this with a probabilistic 
interpretation. Divide all weights Wi and the flow / by the normalizing factor Yl v ev Wl ( v ) 
to obtain new weights Wi and a new flow / which still obeys (1201) and f[ 2 T]) . Now consider 
the following random process. 1 let zq G V be a random variable with probability distribu¬ 
tion given by w\. Given u k G V, I independently choose a u k+1 G N{v k ) with probability 
distribution given by 


P(zq+1 = v\v k ) = f(k, v k , v)/w k (v k ). 
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Notice that for this probability to make sense we would need 

1 = F(u k+1 G N(u k )\u k ) = —]— 'Y f(k,u k ,v), 

w k {yk) . 

veN(v k ) 

which is given by (l20lh Once i\ is chosen for all i G [N ], we have a randomly constructed 
path p G Vn • Let a(p) : = P(p = p) be the probability distribution of p. 

Now observe that for any v G V, by the construction of p and <5, 

wi(v) = P(zq = v) = P(p(l) = v ) = ^ oi(p), 

p(l)=v 

which up to the normalizing factor of w i( v ) establishes (TT6T) for i = 1. For the sake of 

induction, suppose (H6|) is true for some 1 < i < N. Then for any v G V, 

p(^+i = v) = ^2 = u and d+i = v ) 

u£N~ ( v ) 

= ^2 p (d = u)F(u i+ i = v\ui = u ) 

u£N~ ( v ) 

= "22 Wi(u)(f(i,u,v)/wi(u)) 

u^N~ ( v ) 

= ft i ’ U ’ V ) = ™i+l( V ) 

u£N(v) 

by (J2TJ) • But of course we can also write 

P(^i+i = v) = P(p(i + 1) = v) = ^2 & (p)- 

p(i+l)=v 

Thus by induction, and undoing the normalization, we have constructed our desired a. □ 

3. The Skinny Lipschitz Tube Decomposition 

In this section we demonstrate how to combine the bounds from Section |T] and the discrete 
problem from Section [2] to prove Theorem [0 

Proof of Theorem Q Let u t be a solution to the Schrodinger equation with frequency local¬ 
ization (ED- Let /i be given by Proposition [Q H defined as in (TTUT) . and r > 0 given by 
Proposition [2j We will discretize to the lattice Z d x (rZ) C I* 1 x I; for convenience we will 
reserve the letters a, b for general elements of Z d , h will always denote an element of H , and 
x, y will always be used for generic elements of W l . 

We define for n G Z the weight function m : 7L d x (rZ) —y M + by 

m(a, nr) = J \u nT \ 2 p a+H - 
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First observe that m satisfies a conservation law since the translates of /j form a partition of 
unity: 

(22) Y m ( a , nr) = j \u nT \ 2 Y /i a +H = \H\ J \u nT \ 2 n Z d = \H\ J \u nr \ 2 = \H\\\u 0 \\ 2 L 2 . 


a^Jj d 


ae Z d 


Now we would like to show that m also inherits a finite speed property from flFSI) . Let 
A C Z d and let n G Z. We would like to show that 


(23) E m(a,nT) < m(b,(n + 1) 

a£A b£A-\-H 

Expand the LHS using the definition of m: 


T). 


y^m(a, wt) = / \u nT \ 2 y^ji a+H 

aeA J aeA 

/ |^nr| EE H'd+h' 

J A Us- ZJ 


a^A h£H 

A particular a' G A + H may contribute several times in the above integral. Let 

E>i = {a' G A + H;\A n (a' - H)\ > i}, 

so that in other words a' G means there exists at least i distinct choices of h G H such 
that a' — h G A. Then we can write the above sum a different way, obtaining 

\h\ 


J 2m(a,nr) = Y 

aeA t=l ^ 

Now apply the finite speed bound (IFSI) . 


^nr | 


\H\ , 

J2 m (a,nr) <Y I «(■ 

a£A i =1 




Similarly, if we define 


F>i = {b' eA + 2 H- \(A + H) n (b' - H)\ > i}, 


then we can write 


\n\ 

Y m(a, (n + l)r) = Y 

beA+H i= 1 


|^I(n+l)r 



Thus, to prove (|23|) it would suffice to show E>i + H C F>j. Indeed, if 6 G E>i + H then 
6 = a' + h for some a' G A>j and h G if, so by the definition of E>i, 


i < |{(a, h') G A x H] a' = a + h'}\ 

= | {(a, h')eAxH ; b = (a + h) + h'}| 
< |{(c, ti) e(A + H)xH- b = c + h'}\ 
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This last inequality implies b e F> % . We have therefore shown (125]) . 

We are nearing a situation where we can apply Proposition EJ Let R > 0, and let N = 
|~2 R/t] (rounded to the next even number for convenience). We consider the graph on 
G = (Z d , E) on Z d where (a, b) G E if and only if a G b + H. Define the sequence of weights 
Wi : Z^ —y 1R + for 1 < i < N by 

Wi(a) = m(a, (i — N/2)r). 

Now ()23l) implies (TT4l) . and the conservation of mass from (|22l) implies (IT5]) . Thus by Propo¬ 
sition E] there is a weight a : Vn —> M + on the set of paths of length N on G such that 

m(a, (i — N/2 )t ) = S «(p). 

p(i)=a 

Let p be one such path. Choose r = 10d, and let be the piecewise linear path with 
7 p((* — N/2)t) = p{i). By taking V = C^/t, we obtain the velocity bound |V 7 P | < V. 
Consider the tube decomposition function 

f(x,t) = Y Ap) t w(M)- 

p&V N 

Now we verify that f(x, t ) is a pointwise upper bound for \u t (x)\ 2 . Let a G Z d and 1 < i < N. 
Suppose that (x, t) e I^xKis contained in the prism with bottom-left corner (a, (i — N/2)r), 
meaning that 

aj < Xj < aj + 1, (i — N/2 )t < t < (i — N/2 + l)r 

for 1 < j < d. Now let p G Vn be some path with p{i) = a. Then because we chose r 
large enough, and because 7 p ((i — N/2)) = p(i) = a, it follows that (x,t) G T 7p<r . Using this 
observation we can make the bound 

f(x,t) = Y a (p) T i P A x ^) > Y a W = m ( a ’ (* _ N / 2 ) t )- 

PGVn p(i)=a 

Now by (1LCI) and (IF SI) , there exists some absolute constant C such that 
(24) \u t (x) I 2 < C J \u t \ 2 p a < C J |ri(i-iv/ 2 )r| 2 Pa+H = Cm(a,nr). 

Thus for every (x,t) with \t\ < R , 

K 0)| 2 < c Y a (p)A P A x ^)^ 

P&Vn 

which is exactly (with slightly different letters) the bound (I2|) . 

Now we verify ©. Indeed, by (J22 

Y, (10 d) d a(p) < C ®i(o) = C\H\\\u a \\l*. 

P&Vn a£Z d 


□ 




A DECOMPOSITION FOR THE SCHRODINGER EQUATION 


19 


4. Applications to bilinear and multilinear estimates 

In this section we prove the bilinear Strichartz estimate, Theorem [2] and the multilinear 
restriction theorem, Theorem [3] using the tube decomposition. Neither result is new, but I 
think that these proofs illustrate a few interesting points. Perhaps one of the main lessons 
is that neither the bilinear Strichartz nor the multilinear restriction estimates are truly 
dispersive. Indeed observe that the decomposition described in Theorem [1] does not rule out 
the possibility that there is only one tube. Thus by itself it cannot imply any dispersive 
estimates. 

4.1. Proof of Bilinear Strichartz. In this part we prove Theorem [2] using the decom¬ 
position given by Theorem [T] I believe this proof illustrates something intuitive about the 
estimate - roughly speaking the intuition is that a fast and a slow particle cannot be in the 
same place for very long. This intuition is expressed in the proof by observing that a “fast” 
and “slow” tube must have a small region of intersection. 

Recall the definition of the annulus, 

An ■= {£; N/2 < |£| < 21V}. 

We define in addition the annulus 

A* n := {f; N/ 4 < | < 41V}, 

Proof of Theorem^ Let u 0 ,v 0 G L 2 (M d ) have frequency localizations supp u 0 C A N and 
supp vo C Am with M -C N. 

Let V be the constant appearing in Theorem |T] (which is the same constant appearing 
in Corollary [I]), and let ^ C A± be a set of points such that the balls with radius 

1 /(10W) centered at {£*} cover A\ and are contained in A\: 

(100V) d 

Ai C |J B 1 / 10 V (fi) C A* v 

Z— 1 

Using a partition of unity on (a scaled version of) this covering, write the decompositions 

(ioov) d (ioov) d 

Uo = ^2 V 0 = y ^ y j 

i= 1 z=l 

such that supp u t c B N p iou)(Afj) and supp A C P>M/(\QV){Mfj). We can also demand that 
these decompositions are nearly orthogonal in the sense that 

(ioov) d 

\\ui\\li A 2 ||m||^2 
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and likewise with v t . Using this decomposition and the Cauchy-Schwartz inequality, we 
obtain 

(100H) d (100V) d 

IMUf.^xR) < (100VT \\ u i v j\\LlJR d xR)- 

i =1 j=l 

Now we bound these cross-terms lliqi-hll/^ (r^xM)- From Corollary [T] we obtain paths {^ n }Y = i 
with weights {w n }™ =1 such that 

OO 

mm)i 2 <E r ^')n , r N ~ 1 (x,t). 

71—1 

For convenience we will write T n instead of T ln rN -i. We also use Corollary |Tj on v 3 to obtain 
paths { 7 m}m=i with weights {w' rn } < ^ = x such that 

OO 

I Vj(x,t)\ 2 < ^2w' m T! n (x,t) 

772—1 

where T^(x,t) = Ty rM - 1 . Applying these pointwise bounds, we estimate 


(25) / / \ui(x, t)\ 2 \vj(x, t)\ 2 dx dt < I L T n (x, t)T^(x, t ) dx dt. 

Jr JR d m=1 n=1 </R </R d 

Now the paths y n and 7 ^ defining T n and T m have velocities satisfying 


IV 7 n - N£i\ < N /10 


and 

|vy„ - m&i < M/10. 

Since £i,£j G Ai, we can conclude that C A)y and VY m C A)b- Thus the tubes T n 
and T' m can meet at most once, in a region of spatial width N~ x and only for a duration of 
(. MN )“ 1 . That is, 



T n (x,t)T' m (x,t) dxdt < CN~ 1 ~ d M~ 1 . 
Plugging this back into (1251) and then applying (J4j), 

OO OO 


UiV 




772—1 72—1 


< CN^M^^M^w'J^N^Wn) 


772—1 
T —1 71/fd-l IL. II 2 


72—1 


<CN~ X M^ MYMU. 
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Finally applying the almost orthogonality of the decomposition allows us to conclude 

(100V) d (100V) d (100U) d (100V) d 

iMli^xR) < cN - lMd ~\ iKH£g)( 

i —1 j =1 i =1 j =1 

The theorem is finished by taking square roots. □ 

4.2. Multilinear restriction. We can also prove the multilinear restriction estimate using 
the multilinear Kakeya problem for Lipschitz tubes, proved by Guth in {Gut 15] . First let us 
set up some notation. Let {'iy}’L 1 C R" be a collection of unit vectors. We say that {-ty} is 
//-transverse if |ui A V 2 A • • • A v n \ > u. Let T, = be a collection of Lipschitz curves. 

Given <5 > 0, we say that Tj is h-close to parallel with ly if for any tangent vector t of a curve 
7 1 1 — Vj | < 5. Let Tj j denote the tube of unit radius centered along the curve j h j. 

Theorem 5 (Multilinear Kakeya for Lipschitz tubes |Gutl5j b For every e > 0 and v > 0, 
there is a C £ ^ > 0 and a 5 such that the folllowing holds: If {vi} C W 1 is a set of unit vectors 
that is v-transverse, and each T* = {qtj jydj is a set of Lipschitz paths which are 5-close to 
Vi, and Wij is any positive collection of weights, then 

, n / N j \ n / Nt 

( 26 ) / || y / <c^n 

JBr i=\ \j= 1 / i= 1 Vi=l 



E 


Mil?) < 2 CN 


-1 


M' 


d—l | 


U 


12 
I Ll 



This theorem implies a version of multilinear restriction for solutions to the Schrodinger 
equation which are highly localized in frequency. 

Lemma 5. For every £ > 0 there exists a C £ and a 5 > 0 such that the following holds: 
Let r = (10d ) -1 and suppose {^}f =0 is a set of vectors with £ 0 e B r ( 0) and G Brief) for 
1 < i < d. If also Ui are solutions to the Schrodinger equation with supptq G Bs(f,i), then 

p d d 

/ TT x)\ 2 / d dx dt < C e R £ IT \\ui\\% d 

<Rj\t\<Rf = t i=0 

for all R > 0. 

Proof. Consider the vectors v[ = (1,-2^') G M d+1 , and let Vi = v’J |r/'|. The vectors {uj} 
correspond to the unit normals to the paraboloid P = {r = |£| 2 } C M d+1 at the points 
(|£i| 2 > 6 )- Because of the constraint on the positions of G, we can guarantee that the set {iy} 
is //-transverse for some absolute constant v > 0. 

Let e > 0, and let C £ and 5 be given according to Theorem [5] with this constant v and 
with n = dpi. Let V be the same constant that appears in Theorem [IJ and define 5 = 5/V. 
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Now if supp Ui C Bg(^i), then according to Corollary[T] it is possible to find paths 
such that |V 7 ij — £*| <5 and weights {wij} such that 

\ui(t,x)\ 2 < ^2wijT^. trS -i(t,x) 


and 


^(r<5 x )V,i < C\\ui\\ 2 L 2. 


From now on we will simply write T^j for the tube T 7 n 5-1 centered at 7 with width r5 1 . 
Thinking of 7 jj as a curve in M d+1 , where the first coordinate is time and the remaining 
d coordinates are spatial, the condition IV 7 ij — £j| < 5 ensures that the curves { 7 ij} are 
5-close to Vi . Since the radius of the tubes T t J is rd -1 , which is much larger than 1, we may 
by splitting the tubes into thinner pieces think of them as having radius 1 (buying simplicity 
at the cost of an ignorable factor of 5 £ ). Thus we may apply Theorem [5] to obtain 

1/d 


'\x\<R/2 J\t\<R/2 i=i 


d p d / 00 \ 

JJ \Ui(t,x)\ 2/d < JJ I J2 w i,i T hj ) 
i=0 J B R i =0 \j =1 / 


00 \ 

< c,r- n e w >j 

i=0 \j=l 


< 


c , r ‘ niKn« 


»=0 


□ 


The multilinear restriction theorem, Theorem [3j follows upon decomposing each of the 
d + 1 solutions into finitely many pieces with frequency localizations given by 5, and then 
summing the finitely many contributions. 


References 

[BCT06] J Bennett, A Carbery, and T Tao, On the multilinear restriction and Kakeya conjectures , Acta 
mathematica 196 (2006), no. 2, 261-302. 

[BD14] J Bourgain and C Demeter, The proof of the l 2 decoupling conjecture, arXiv preprint 
arXiv:1403.5335 (2014). 

[BG11] J Bourgain and L Guth, Bounds on oscillatory integral operators based on multilinear estimates, 
Geometric and Functional Analysis 21 (2011), no. 6, 1239-1295. 

[Bou98] J. Bourgain, Refinements of Strichartz’ inequality and applications to 2D-NLS with critical non¬ 
linearity, Intern. Mat. Res. Notices 5 (1998), 253-283. 

[CKS + 02] J. Colliander, M. Keel, G. Staffilani, H. Takaoka, and T. Tao, Almost conservations laws and global 
rough solutions to a nonlinear Schrodinger equation, Math. Res. Letters 9 (2002), 659-682. 

[FF56] LR Ford and DR Fulkerson, Maximal flow through a network, Canadian journal of Mathematics 
8 (1956), no. 3, 399-404. 


A DECOMPOSITION FOR THE SCHRODINGER EQUATION 


23 


[Gut 15] 
[Hanl2a] 
[Hanl2b] 
[Jos07] 


L Guth, A short proof of the multilinear Kakeya inequality, Mathematical Proceedings of the 
Cambridge Philosophical Society, vol. 158, Cambridge Univ Press, 2015, pp. 147- 153. 

Z. Hani, A bilinear oscillatory integral estimate and bilinear refinements to Strichartz estimates 
on closed manifolds, Analysis and PDE 5 (2012), 339-362. 

Z. Hani, Global well-posedness of the cubic nonlinear Schrodinger equation on closed manifolds, 
Communications in Partial Differential Equations 37 (2012), no. 7, 1186-1236. 

S. Joseph, The max-flow min-cut theorem, http: //www .math .uri . edu/~eaton/maxf lowmincut .pdf. 


2007. 


[KRT02] S. Klainerman, I. Rodnianski, and T. Tao, A physical space approach to wave equation bilinear 
estimates, Journal dAnalyse Matliematique 87 (2002), no. 1, 299-336. 

[SM] A Staples-Moore, Network flows and the max-flow min-cut theorem, 
http://www.math.uchicago.edu/~may/VIGRE/VIGRE2009/REUPapers/Staples-Moore.pdf 
[TaolO] T. Tao, A physical space proof of the bilinear Strichartz and local smoothing estimate for the 
Schrodinger equation, http://terrytao.files.wordpress.com/2013/07/bilinear.pdf, 2010. 
E-mail address: felipeh@mit.edu 


